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q , PACS. 74.20Rp - Pairing symmetries (other than s-wave). 

Q ■ PACS. 74. 25. Ha - Magnetic properties. 

PACS. 74.25.Fy - Transport properties (electric and thermal conductivity, thermoelectric ef- 
O ,■ fects, etc.). 
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Abstract. - Broken time reversal symmetry (BTRS) in d wave superconductors is studied and 
is shown to yield current carrying surface states. The corresponding spontaneous magnetization 
5> is temperature independent near the critical temperature T c for weak BTRS, in accord with 
recent data. For strong BTRS and thin films we expect a temperature dependent $ with a 
paramagnetic anomaly near Tc- The Hall conductance is found to vanish at zero wavevector q 
ii 1 and finite frequency u), however at finite q,oj it has an unusual structure. 

(N 
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£■ — ' Recent data on the high T c superconductor YBaiCu-sO x (YBCO) has supported the pres- 

ence of broken time reversal symmetry (BTRS) 0-0] . A sensitive probe of BTRS are Andreev 
surface states. For a d wave with time reversal symmetry bound states at zero energy are ex- 
pected for a surface parallel to the nodes (a (110) surface in YBCO). When BTRS is present, 
by either a complex order parameter or by an external magnetic field, the bound states shift 
to a finite energy. Indeed tunneling data usually shows a zero bias peak which splits in an 
applied field; the splitting is nonlinear in the magnetic field, indicating a proximity to a BTRS 
state JMW]. In fact, in some samples tunnelling data shows a splitting even without an external 
field |lj,|2|, consistent with BTRS; the splitting increases with increasing overdoping (|,||], 

Further support for a spontaneous BTRS state are spontaneous magnetization data as 
observed in YBCO |l, setting in abruptly at T c and being almost temperature (T) independent 
below T c . The phenomenon has been attributed to either a d x 2_ y 2 +id xy state or to formation 
of 7r junctions. No microscopic reason was given, however, for the spontaneous magnetization 
being independent of both T and of film thickness B . 

It has been shown theoretically that BTRS can occur locally in a d x i_ y i superconductor 
near certain surfaces P-pTfl leading to surface currents. The onset of such BTRS is expected 
to be below T c and therefore does not correspond to the spontaneous magnetization data H. 
We note that in response to an external magnetic field the surface states are paramagnetic 
and compete with Meissner currents. This effect has been proposed to account for a minimum 
in the magnetic penetration length p]. In fact, it was proposed that this paramagnetic effect 
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leads to a spontaneous currents and BTRS in a pure d x 2_ y 2 state J7lp[. The onset of this 
BTRS is much below T c M and therefore does not correspond to the data 0] . 

Of further theoretical interest is the relation of the BTRS state to quantum Hall systems 
with a variety of Hall effects [12] [l5[ . In particular a finite charge hall conductance has been 
suggested |12[, though this has been questioned [ pT[ . 

In the present work we assume that BTRS is a bulk property, i.e. both components of 
an order parameter d x 2_ y 2 + id xy set in at T c ; this is a plausible scenario which can possibly 
account for the data. In previous works |9J-]ll[| surface states appear due to a surface induced 
id xy component. We find, however, the less anticipated result that the bulk state d x 2_ y 2 +id xy 
leads to surface states with finite surface current densities. The latter situation was found in 
the bulk p wave state [[16| and for the total current was inferred from topological considerations 
P7|| . We then evaluate the spontaneous magnetization and show that it is dominated by (100) 
surfaces; for thin films it increases with the ratio A'/A (A and A' are the amplitudes of d x 2_ y 2 
and d xy , respectively) while for thick films it has a maximum at A/£' ~ 1 where £, £' are the 
coherence lengths associated with A, A', respectively (£' =Fermi velocity/ A', similarly for £) 
and A is the magnetic penetration length. The maximum for YBCO is at A'/A rs £/A rs 0.01. 
Throughout we assume an exterme type II superconductor, i.e. £ <C A, while £'/A is arbitrary. 
We show that for weak BTRS, A/£' < 1, the spontaneous magnetization is T and thickness 
independent, while for strong BTRS thickness and T dependence may occur, as well as a 
transition to a paramagnetic state close to T c . For the sample of Ref M we estimate A'/A « 
10 -4 , i.e. weak BTRS. We also derive the effective action and identify the Hall coefficient 
which has an unusual wavevector and frequency dependence. 

First we demonstrate the existence of surface states. Consider a d x 2_ y 2 + id xy state where 
the order parameter is 

A(p x ,p v ) = A'p x p y /kp + iA(f x - f v )/k F (1) 

where p = — i?iV is the momentum operator and kp is the Fermi momentum. We consider 
a vacuum-superconductor boundary at x=0, and assume for now that A, A' are constants at 
x > and vanish at x < 0. For A ^> A' this corresponds to a (100) surface; to describe a (110) 
surface A and A' need to be interchanged. The electron-hole wavefunctions u(x, k y ) exp(ik y y) , 
v(x, k y ) exp(ik y y) with mass m satisfy the Bogoliubov de-Gennes equations 

(2m) - (—k F — d jdx )u{x,k y ) + A(p x ,k y )v(x,k y ) = eu(x,k y ) 
A*(p x ,k y )u{x,k y ) + (2my 1 (kp + d 2 /dx 2 )v(x,k y ) = ev(x,k y ) (2) 

The decaying eigenfunctions have the form ~ exp(±ifca; — x\k y \/t;'kp) with k = \k\ and 
£' = kp/mA'. Specular reflection at the surface requires a superposition of ±fc states which 
vanish at the surface. This yields the eigenvalue equation 



y/\A(+k,ky)\*-e* _ A(+k,k y ) 



-ie+^\A(-k,k y )\ 2 -e 2 A(-k,k y ) 

Its solutions are e = —sign(k y )A(k 2 — k y )/k F . In terms of the incidence angle £, k y = kp sin (, 
k = kp cos C, the allowed positive eigenvalues are e = A cos(2£) for — 7r/4 < £ < and e = 
—A cos(2£) for 7r/4 < C < 7r /2 (see inset of Fig. 1). We note that self consistency would imply 
that A' = at x = |l|]; the eigenfunctions would then be ~ exp[— J Q X A'(x')dx'\ sin(\/vp], 
resulting in a qualitatively similar dependence on £'. The dominant order parameter A is 
finite at the (100) boundary [pTj, hence we expect our results to be quantitatively correct. 
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Fig. 1 - Spontaneous flux for a (100) boundary in thick films (£' < d) in units of 2(f>oL y \£±./-K\^ ) T c . 
The hatched areas in the inset show the allowed directions (k,k y ) of bound states. Their spectra 
spans the range e = (diagonal lines) up to e = A as shown. 



In presence of a vector potential A y (x) the spectrum e is Doppler shifted by ~(e/mc)k y A y . 
For the expectation value of the current density in the y direction and the charge density we 
obtain 



Air 

Jedge\%) 



Tledge \p^J 



-% f /2 dC cos C sin 2 Ce-W/? tanh( Ac ° s2C ± (*/<** ^W*)) 
n £ ^o Jo 2T 



ekp 
^d?Jo 



tt/2 



<i£cos£sin£e' 



-2a:sinC/5' 



(1) 



where Ao = (mc 2 d/2k 2 7 e 2 ) 1 ^ 2 — X(T = 0) and d is the spacing between layers, converting the 
current per layer of the states in Eq. (pi) to a current density. Note that for either A = 
or A' — all angles £ are allowed in the solution of Eq. (ph and then the current density 
vanishes. This demonstrates then that BTRS leads to current carrying surface states. We 
note also that the integrated current J °° j e d ge (x)dx vanishes, unlike the p wave case fig] . 

The response to j e dge(x) includes London terms as well as a BTRS induced Chern-Simon 
term. As shown below, the Chern-Simon term is weaker at T — > T c so that London's equation 
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with jedge{%) as a source term is valid, 

- V 2 A y {x) = [-(l/X 2 )A y (x) + (An/c)j edge {x)]e{x) (5) 

where 9{x) is a step function. We consider first solutions where the Doppler shift ~ A y on 
the right hand side of Eq (Q) is neglected. This is solved by the Greens' function G(x,x') = 
— (A/2)[cxp(— \x — x'\/X) + exp(— | a; + x'|/A)] resulting at T — ► T c in the solution 



,tt/2 

(0) = (20oAA/7rAoT c ) / dCcosCsin 2 Ccos2C(2sinC + £'/A) _1 (6) 

Jo 



The total spontaneous flux is <J> = A y (Q)L y where L y is the length of the boundary. The ratio 
§ = -$/(2(j>oLy\A/Tr\lT c ) is shown in Fig. 1. It varies between £'/12A at A > £' to A/15£' 
at A < £' with a maximum of 0.014 at A « £'. For a (110) surface l> = £A'/12AA, much 
smaller than for a (100) surface. We note that $ is T independent up to T c since the product 
A A is finite at T — ► 0, consistent with the spontaneous magnetization data |g|; more details 
on the data follow in the discussion. 

The above refers to thick films with thickness d > £'. For thin films the right hand side of 
Eq. (j5|) is multiplied by dS(z). For A y (a;) = A y (a;, z = 0) this yields 

4,(3) - A y (0) - / da/ In | (x - a/)/a/|[-l/A 2 )Ay(0 + (^/c)j edge (x')]d (7) 

Jo 

which implies a slow (non-exponential) decay of A y (x). To avoid divergence of dA y (x) / dx\ x =o 
the relation A a (0) = A 2 (4ir / c)j e d ge (0) must hold. This, interestingly, yields for $ the previous 
result of the A <C £' case, i.e. $ = A/15£'; in Fig. 1 this is the tangent line to the thick 
film curve at the origin. Hence we can define two regimes: Weak BTRS with A/£' < 1 where 
the spontaneous flux is T and d independent, and strong BTRS with A/£' > 1 where film 
thickness matters. In the latter case T dependence is induced as £' < d changes to the thin 
film case £' > d as T — > T c . 

We now reconsider the effect of the Doppler shift. Near T c the effect is linear in A y and 
can be incorporated as an effective A e // replacing A in the results above, where 

J- ~ 1 _ £° (8 ) 

X 2 eff ~ A 2 A 2 max(e,A) ' W 

Hence at a temperature T^ A e // changes sign, where (T c — T' c )/T c rs [^/max(f, A)] 2 , i.e. less 
then 10 -4 for YBCO. The effect of the Doppler shift was previously considered only at low 
temperatures where on the (110) surface paramagnetism and spontaneous surface currents 
appear even when A' = at ~ {£,o/Xq)T c <C T c |7|,|j. For thin films max(£',A) is replaced 
by £' and (T c - T' c )/T c w (£/£') 2 also for strong BTRS where T c - T' c is enhanced. Hence a 
transition from paramagnetic (T > T' c ) to diamagnetic response (T < T' c ) can be observed at 
T' c and the magnetization <l = (A/15£')(A e ///A) 2 is peaked and changes sign at T' c . The effect 
is more pronounced for the (110) surface, however, the condition for a thin film d < £ is more 
difficult to achieve. 

Next we consider the effective action of a bulk d x i_ y i + id xy superconductor. In terms 
of the Nambu spinors ^ ( r ) = [u*(r),«*(r)], the superconducting phase 0(r) and Pauli ma- 
trices t,, the transformation ip(r) — > exp[iT30(r)/2]y>(r) yields the off-diagonal Hamiltonian 
/ d 2 rip^ (r)/iA'0( r ) where 

ft A = -[A(-9 2 + 9 2 )n + A'a^raj/A; 2 (9) 
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and we neglect terms with V# << kp. The action in presence of electromagnetic potentials 
A, tp is then 

S = f d 2 rdt^ (id t - r 3 e(p) - h A - £)V 

S = r 3 (a + a 2 /2m) + a ■ p/m - iV ■ a/2m (10) 

where e(p) = (p 2 — k F )/2m and we introduce the gauge invariant potentials a = ^V# — eA 
and a — \-§^0 ~ eip. Expansion to 2nd order in a, ao leads to the effective action 

/j2 j 
— ^F^(q,w)a^(q,cj))a^(-q,w)) (11) 

At T = and q, to — ► we obtain Poo = No (density of states which is JVo = to/27t in two 
dimensions), Py — — A^c 2 where c 2 s = vf/V%, while Poj(q) = isign(AA')eoijqi/(4:n) and eoij 
is the antisymmetric unit tensor. The latter term reflects BTRS and is derived for A' <C A. 

Integrating out the phase we obtain the effective action in terms of the electromagnetic 
potentials A, ip 

/d oduj co i 

7^p-{ ^(J^ t-PooMq^)! 2 " ^eoy%v(q,w)^j(-q,-w) 

+0(c 2 |A| 2 )]-Poo(^) 2 |A(q, w )| 2 } (12) 

The total electromagnetic action includes also the Maxwell field part Sm — J d 2 rdt(E 2 — 
H 2 )/8ir. E.g., for u> ^ 0, q — > the propagator for ip yields the plasmon mode at ui p — c/Xq. 
The Hall current J y is identified by a functional derivative with respect to A x leading to the 
Hall coefficient 

2 2 2 

a xy (d,uj) = sign(AA') e C f (13) 

Transport is defined by taking the q — > limit first, i.e. a xy — 0. Hence the conventional 
Hall coefficient vanishes, as expected from Galilean invariance |Q. A limit in which u — > is 
taken first yields a quantized "static" conductance e 2 /2h which was argued to correspond to 
o'xy 7^ in presence of a boundary |0|. In absence of an external magnetic field, and given 
a spontaneous magnetization decaying in the bulk (as confirmed below), Ampere's law yields 
zero total current, hence a xy — 0; this is valid also with a boundary and external electric 
field. It is intriguing, however, that a xy {(\,oj) has a nontrivial structure and space resolved 
measurement of a Hall current could then probe the full Eq. (|l|). We note that a result 
similar to Eq. (|l3| ) was obtained for superfluid 3 Hc Jig] , 

We proceed to derive the effective action in presence of a boundary and at finite T. Spe- 
cial care is needed for the Chern-Simon coefficient Pqj which now nonlocal due to specular 
reflection at the boundary. After integrating out 9 we obtain the form (considering only the 
u) = term) 



. c s 



dA y (v) 



S,,{A.-p\ -- e 2 J dr{P QQ (p 2 (r) - (^) 2 A 2 (r)) +b 1 (rMr)A y (r) + b 2 (rMr)^j^}il4' ! 



c 



Variation of S + Sm leads to a generalization of Eq. (H) in which the equations for A and ip 
are coupled, 



1 d 2 , 47re 2 , dA y 4ne 2 



( a! " b* )v ~ ~^ h2 ^ ~ -rfT blAy = 4 ^^^ ( 15 ) 



VJ 
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, 1 d 2 . . Aire 2 , <9w Aire 2 ,, „ , . 47r . . , „. 

( A 2 " " &?^ " ^T 52 fe + ~W {bl - Wlf = -T 3ed3e{x) (16) 

where A^ = (87re 2 iVo) _1 / 2 is the Debye screening length. For T — + T c we obtain A « Ao(l — 
T/Tc)- 1 / 2 , bi(r) = Q.ll(AA'/T?)£ln[AA']/2hcd and 6 2 (r) = 0.21(AA'/T 2 )/2W. 

An external electric field leads to ip ~ exp(— x/\<j), hence the magnetization remains 
localized even in presence of such a field and space integration of Eq. Jig ) (i.e. Ampere's 
law) yields a zero total current, i.e. a xy = 0. The Chern-Simon term, however, affects the 
spontaneous magnetization, leading to an additional flux ~ (A'/A)(l — T/T c ) which vanishes 
at T -> T c . 

We consider now in more detail the experimental data on the spontaneous magnetiza- 
tion ||. The data shows that for a YBCO disc with a perimeter of L y « 2cm the spontaneous 
magnetization is temperature independent in the range 80-89K and is also thickness indepen- 
dent in the range 30-300nm with a value of ~370o- Taking A A « AoAo, their T — value, 
and typical YBCO parameters we find $ rs 10~ 3 . Fig. i implies that the limit £' > A applies 
i.e. weak BTRS; hence for either thick or thin films we estimate A/£' sw 10~ 2 or A'/A « 10~ 4 . 
We propose therefore that increasing the ratio A'/A, e.g. by using overdoped YBCO H, one 
can enhance the spontaneous magnetization up to a maximum of sw f 3 </>o when A'/A w O.Of 
within the thick film regime. 

For strong BTRS, A/£' > f, the film thickness matters, i.e. we expect a temperature 
dependence due to the crossover from thick to thin film regimes at d « £' as T — + T c . For thin 
films (d <£ < A) we obtain l> = A/12£', i.e. for YBCO the total flux can reach 10 5 A'/A</> 
per cm of boundary, much higher than thick film values. The situation of a strong BTRS with 
thin films is interesting also as being the most likely one to show the paramagnetic anomaly 

atT^T c [l-(e/e') 2 ]- 

In conclusion, we have shown that surface states of a d x i_ y i + id xy superconductor lead 

to spontaneous magnetization which is T independent and thickness independent for weak 

BTRS, X/C < I, in accord with the data ||. For strong BTRS, A/£' > I, as expected in 

overdoped YBCO , a crossover from thick to thin film behavior can lead to T and thickness 

dependence, as well as to an observable paramagnetic anomaly near T c . We also find that 

the Hall conductance has an unusual a xy (q, w) dependence, though its conventional transport 

value vanishes. 
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